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We report the first observation of J/ip riK*^K*'^ decay in a J/ip sample of 58 million events collected with 
the BESII detector. The branching fraction is determined to be (1.15 ± 0.13 ± 0.22) x 10^^. The selected signal 
event sample is further used to search for the F(2175) resonance through J/i/> ^ 77r(2175),F(2175) ^ K*°K*°. 
No evidence of a signal is seen. An upper limit of Br( J/V^ ?7F(2175)) • Br(r(2175) ^ /s:*°¥*'^) < 2.52 x 10"" 
is set at the 90% confidence level. 



1. Introduction 

Following the observation of ^(2175) by the 
BaBar Collaboration in e+e~ 77SK0/o(98O) 
via initial-state radiation ilj, the resonance was 
observed by the BES Collaboration in J/ip 
r]4>f o{980) [2] and more recently by the Belle Col- 
laboration in e+e~ ^ 77s_R07r+7r~ [3]. Since both 
the y(2175) and y(4260) [4] are observed in 
e~^e~ annihilation via initial-state radiation and 
these two resonances have similar decay modes, it 
was speculated that y(2175) may be an s-quark 
version of y(4260) [1]. There have been a num- 
ber of different interpretations proposed for the 
y(4260), that include a gcc hybrid ^5] [3] [7j, a 
A^Si cc state [S], a [cs]5'[cs]s tetraquark state 
[9], or a baryonium [10]. Likewise y(2175) has 
been correspondingly interpreted as: a gss hybrid 
[IT] , a 2^Di ss state [12], or a ssss tetraquark 
state [13]. None of these interpretations has ei- 
ther been established or ruled out by experimen- 
tal observations. 

According to Rcf. [12], a hybrid state may 
have very different decay patterns compared to a 
quarkonium state. Measuring the branching frac- 
tions of some decay modes may shed light on un- 
derstanding the nature of F(2175). Among those 
promising decay modes, r(2175) K*°K*^ is of 
special importance. This decay mode is forbidden 
if y(2175) is a hybrid state but allowed if it is a 
quarkonium state. 

On the other hand, there are still lots of 
unknown decay modes of J/ij) and investigat- 
ing more of them is useful to understand the 
mechanism of J/V' decays. Based on a sam- 
ple of 58M J/V' events collected by the BESII 
detector at the Beijing Electron-Positron Col- 
lider (BEPC), a search for the process J/i) — > 
77^(2175), r(2175) K*°'K*° is performed. In 
addition, the first measurement of the branching 
fraction Br( J/V' riK*'^K*'^) is obtained. 

2. Detector and data samples 

The upgraded Beijing Spectrometer detector 
(BESII) was located at the Beijing Electron- 
Positron CoUidcr (BEPC). BESII was a large 
solid-angle magnetic spectrometer which is de- 



scribed in detail in Ref. [14]. The momentum 
of charged particles is determined by a 40-layer 
cylindrical main drift chamber (MDC) which has 
a momentum resolution of Cp/p=1.78%'\/l -\- 
{p in GeV/c). Particle identification is accom- 
plished using specific ionization (dE/dx) mea- 
surements in the drift chamber and time-of-flight 
(TOF) information in a barrel-like array of 48 
scintillation counters. The dE/dx resolution is 
^dEjdx — 8.0%; the TOF resolution for Bhabha 
events is (Jtof = 180 ps. Radially outside of 
the time-of-flight counters is a 12-radiation-lcngth 
barrel shower counter (BSC) comprised of gas 
tubes interleaved with lead sheets. The BSC mea- 
sures the energy and direction of photons with 
resolutions of as/E ~ 21%/^/E {E in GeV), 
(70 = 7.9 mrad, and Uz = 2.3 cm. The iron flux 
return of the magnet is instrumented with three 
double layers of proportional counters that are 
used to identify muons. 

A GEANT3 based Monte Carlo (MC) pack- 
age (SIMBES) [T5] with detailed consideration of 
real detector performance (such as dead electronic 
channels) is used. The consistency between data 
and Monte Carlo has been carefully checked in 
many high purity physics channels, and the agree- 
ment is quite reasonable |15j . 

3. Analysis 

The decay channel under investigation, J/ij} 
riK*^'K*°, r] 77, K*° if+Tr", ¥*° K~Tr+, 
has two charged kaons, two charged pions, and 
two photons in its final state. A candidate event 
is therefore required to have four good charged 
tracks reconstructed in the MDC with net charge 
zero and at least two isolated photons in the BSC. 
A good charged track is required to (1) be well 
fitted to a three dimensional helix in order to en- 
sure a correct error matrix in the kinematic fit; 
(2) originate from the interaction region, i.e. the 
point of closest approach of the track to the beam 
axis is within 2 cm of the beam axis and within 20 
cm from the center of the interaction region along 
the beam line; (3) have a polar angle 9, within the 
range | cos6'| < 0.8; and (4) have a transverse mo- 
mentum greater than 70 MeV/c. The TOF and 
dE I Ax information is combined to form a particle 
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identification confidence level for the tt, K, and p 
hypotheses, and the particle type with the high- 
est confidence level is assigned to each track. The 
four charged tracks selected are further required 
to be consistent with an unambiguously identi- 
fied K^TT^ K^TT^ combination. An isolated neu- 
tral cluster is considered as a good photon when 
(1) the energy deposited in the BSC is greater 
than 60 MeV, (2) the angle between the near- 
est charged track and the cluster is greater than 
15°, (3) the angle between the cluster develop- 
ment direction in the BSC and the photon emis- 
sion direction is less than 30° , and (4) at least two 
layers have deposits in the BSC and the first hit 
is in the beginning six layers. A four-constraint 
(4-C) kinematic fit is performed to the hypothe- 
sis J/i/j 'yjK~^K~7r~^Tr~ , and if there are more 
than two good photons, the combination with the 



of good photons to be two. 



smallest x: 



value is selected. We fur- 



ther require that X-,-,^+jf-^+^- ^ 20. Because 
we are not interest in the events of which the 
two photons come from ir^, we require the invari- 
ant mass of two photons to be greater than 0.3 
GeV/c2. 

3.1. Branching fraction of J/i) r]K*°'K*° 

After applying the above event selection cri- 



teria. Fig. 1(a) shows the scatter plot of M/f+^- 
versus Mk-^+. One can see K*°'K*°, K*°K-Tr+, 
K*'^K~^n~ , and K'^n" K~tt'^ events scattered 
in different regions of the plot. The signal re- 
gion in this analysis is defined by |M^±jrT — 
TO;f.o(m^,o)| < 0.05 GeV/c^, which is shown as 



the middle box in Fig. 1(a) Other boxes shown 



are side-band regions, and events in these regions 
are used to estimate the background in the sig- 
nal region. The K^n^ invariant mass spectra are 



shown in Fig. 1(b) where the solid histogram is 
K^TT^ and the dashed histogram is K'tt'^. Fig- 
ure 2(a) shows the 77 invariant mass spectrum for 



events in the signal region, where an 77 is seen. In 



Fig. 2(a) the shaded histogram is the spectrum 



obtained requiring two good photons, while the 
dashed histogram is the spectrum for more than 
two photons. When there are more than two pho- 
tons, the ratio of signal over background is much 
lower. In order to remove potential backgrounds 
as much as possible, we also require the number 



Figure 2(b) shows the 77 invariant mass spec- 



trum of events surviving the above selection, 
while the shaded histogram is the normalized 
background estimated using the side-band regions 
shown in Fig. 1(a) The number of J/rp 
r]K*^K*'^ events is determined by fitting the spec- 
tra in Fig. [2(b)| The J/'ip ■qK*^'K*^ branching 



fraction is determined using 



■e-Br{K'-0->K+TV-)-Br{K*O^K-TT+ )-Br{T]^-f-f) ' 



where Ngig = 347 is the number of events in the 
signal region, obtained by fitting the spectrum 
in Fig. 2(b) (the blank histogram); Nsb = 138 is 



the number of background events estimated from 
side-band regions, obtained by fitting the spec- 



trum in Fig. 2(b) (the shaded histogram); iVj/^ is 
the total number of J/ip events [T7]; e = 1.79% is 
the detection efficiency obtained from MC simula- 
tion_of J/i; T]K*"'K*°; and Br(7f *° ^ K+n-), 
Br{K*'^ K'TT^) and Br{ri 77) are the cor- 
responding branching fractions. Figures 3(a) and 
|3(b)| show respectively the fitting results of the 
signal and side-band events, where the shape of 
the 77 invariant mass spectrum obtained from 
the MC sample J/^ riK*"^*^ is used as the 
signal shape and a third order Chebyshev poly- 
nomial is used as the background shape. The 
J/^p riK*^K*^ branching fraction is deter- 
mined to be 



B 



'{J/jP riK'^K*") = (1.15 ± 0.13) x 10-^ 



where the error is statistical only. It is the first 
measurement for this decay mode of J/ip and it 
is shown that this mode is a typical three bod- 
ies decay. The branching fraction is compatible 
with the result of Br{J/ip r]K~^ K~Tr'^Tr~) = 
(1.84±0.28) X 10-3 given by BaBar Collaboration 
[16] . It is worth mention of that this branching 
fraction is several times smaller than the radia- 
tive decay mode J/ip jK*^K*'^ which is very 
different from the situation of pp that the branch- 
ing fraction of J/ip — > r]pp is much bigger than 
J/ip IPP- 
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Figure 1. (a) Scatter plot of Af^+^- versus M;f-^+ invariant mass, where the middle box is the signal 
region and the other boxes are the side-band regions, (b) The invariant mass spectra of K^tt^; the solid 
histogram is K~^Tr~ and the dashed is K~7t'^. 



3.2. J/tp 77^(2175) ^ r]K*°K*° 

Next, we search for a possible resonance re- 
coiling against 77. So in addition to the above 
requirements, we require that the 77 invariant 
mass satisfies — rn^| < 0.04 GeV/c^ and 

define the side-band region to be O.f GeV/c^ 
< \M^^ - mn\ < 0.14 GeV/c^. The K*°K*° 
invariant mass spectrum recoiling against 77 for 
J/ip j]K*°'K*'^ is shown in Fig. gl where the 
dashed histogram is the contribution from phase 
space for J/V' riK*^K*'^ and the shaded his- 
togram is the contribution from the normalized 
side-band events in the 77, K*'^ and K*'^ side-band 
regions. There is no obvious enhancement in the 
region around 2.175 GeV/c^. 

The backgrounds in the selected event sample 
are studied with MC simulations. For the decay 
J/i(j r]K*'^K*'^, the possible main background 
channels are: J/iJj -> r/K*°'K*° {3tt°)K*°'K*°; 
J/^Jj a+K-K*^ {r]n+)K-K*° + c.c; J/^) 
p+K*-I£° {tt+tt°){K-7t°)K*° + c.c; J/ijj 
'yTr°K*°'K*°; J/ij; (j)r]' K+ K-'qir+TT- ; for 



each channel a sizable MG sample is simulated. 
There is no peak around 2.175 GeV/c^ in the 
K*^K*^ invariant mass distribution in any back- 
ground channel. 

We fit the mass distribution to determine a 
possible signal, where three parts are included in 
the total probability distribution function (p.d.f): 
(1) for the signal p.d.f, we use the shape of the 
j{*Oj^*o invariant mass spectrum obtained from 
MC simulation of J/ip r]Y{2175) -> riK*^'K*^ 
produced with the mass and width of y(2175) 
fixed to BaBar's results; (2) for the normalized 
phase space contribution p.d.f., we use the shape 
of the if*°_ftr*° invariant mass distribution ob- 
tained in the J/-0 t]K*°'K*° MG simulation, 
normalized with the branching ratio obtained in 
the previous section; (3) for the other possible 
backgrounds, we use a third order Chebyshev 
polynomial. 

The product branching ratio is determined us- 
ing 

Br{J/Tp riY{2175)) ■ Br(F(2175) K*°'K*°) = 
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iVjy^ ■e-Br(Jf *°->if+7r-)-Br(7f •0^if-7r+) Sr{77^77) 

(0.7 ±0.8) X 10"*, 

where N°'"' — 11 ± 12 is the number of sig- 
nal events, Nj/^ is the total number of J/ip 
events |T7], e = 1.57% is the detection effi- 
ciency obtained from MC simulation of J/ip — s- 
77^(2175) r]K*°'K*", where the first step de- 
cay used an angular distribution 1 + cos^ 9, 9 
is the polar angle of the rj momentum in the 
center of mass frame, Br{K*^ K~^Tr~) and 
Br{K*'^ if^7r+) and Br{ri — > 77) are the 
corresponding branching fractions. The error is 
only the statistical error. The signal significance 
is only 0.88(7. 

The upper limit of Br{J/Tp r;r(2175)) • 
Br(y(2175) -> K*°K*°) at the 90% confidence 
level is obtained using a Bayesian approach |18j . 
We obtain the upper limit: 

Br{J/iP 77^(2175)) ■ Br{Y{2l75) K''°'K*°) < 


N,/^-€-Br{K'<^^K+iT-)-Br{K*O^K-w+)-Br{r]^-y-f)-{l-a^y^) 

2.52 X 10-", 

where N°^^ = 31 is upper limit at the 90% con- 
fidence level, cr*^'* is the systematic error dis- 
cussed below, and the other symbols are defined 
as above. 

4. Systematic Errors 

In this analysis, the systematic errors on the 
branching fraction and upper limit mainly come 
from the following sources: 

4.1. MDC Tracking efficiency and kine- 
matic fitting 

The systematic errors from MDC tracking and 
kinematic fitting are estimated by using simula- 
tions with different MDC wire resolutions [15] . 
In this analysis, the systematic errors from this 
source are 12.8% for J/?/' ^ riJ{*°'K*° and 12.0% 
for J/V- 77^(2175) r]K*°K*°. 

4.2. Photon detection efficiency 

The photon detection efficiency is studied in 
reference [TS]. The results indicate that the sys- 
tematic error is less than 1% for each photon. 
Two good photons are required in this analysis, 



so 2% is taken as the systematic error for the pho- 
ton detection efficiency. 

4.3. Particle identification (PID) 

In references [TS] and [12], the efficiencies of 
pion and kaon identification are analyzed. The 
systematic error from PID is about 1% for each 
charged track. In this analysis, four charged 
tracks are required, so 4% is taken as the sys- 
tematic error from PID. 

4.4. Uncertainty of intermediate decay 

The branching fraction uncertainties for rj 
-f-f and K*°(K*°) ^ K+tt-{K-tt+) from PDG08 
[18] are taken as systematic errors. 

4.5. Number of J/i/j events 

The number of events is (57.70±2.62) x 10^ 
determined from the number of inclusive 4-prong 
hadrons [T7]. The uncertainty 4.72% is taken as 
a systematic error. 

4.6. Fitting 

4.6.1. J/ijj r]K*°'K*° branching fraction 

When fitting the 77 invariant mass spectrum, 
as described in section III. A, the r] signal shape 
obtained from MC is fixed, and different order 
polynomials are used for the background shape. 
The difference is taken as the systematic error for 
the background uncertainty. We also use differ- 
ent regions in fitting the invariant mass spectrum. 
The total systematic error from fitting is 6.7%. 

4.6.2. Br(J/V' r;r(2175)) • Br(r(2175) ^ 
K*°K*°) upper limit 

When fitting the invariant mass spectrum of 
K*^K*°, as described in section III.B, there are 
three sources of systematic error: for the first 
p.d.f, we used the different resonance parameters 
measured by BaBar and BES, and take the differ- 
ence as the systematic error from the uncertainty 
of signal parameters; for the second, the system- 
atic error comes from the error of the branching 
fraction of J/ip rjK*^K*'^ measured in section 
III. A; for the third, we used the difference be- 
tween fitting with a third order Chebyshev poly- 
nomial and fitting with the invariant mass shape 
from K*^K*'^ side-band events as the systematic 
error for the background uncertainty. Combin- 
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Figure 2. (a) The 77 invariant mass spectrum for data; the dashed histogram is from the > 2 events, 
the shaded histogram is from the — 2 events, and the blank histogram is from all events, (b) The 
77 invariant mass spectrum for Nj = 2, where the blank histogram is from signal region events, and the 
shaded one is from the side-band regions events. 
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Figure 3. Unbinncd fitting results of 77 invariant mass spectra: (a) for the signal region events; (b) for the 
side-band region events, where the signal shape is obtained from the MC 77 invariant mass distribution 
and the background shape is a third order Chebyshev polynomial. 




Figure 4. The K*^K invariant mass spectrum, where points with error bars are candidate events, the 
dashed histogram is from MC phase space for J/-0 riK*^K*'^, the shaded histogram is from side-band 
events, and the solid curve is the fitting result, where the Y(2175) shape used is from MC simulation. 
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ing these contributions, 16.3% is obtained as the 
systematic error from fitting. 

4.7. Different selection of side-band re- 
gions 

We used different side-band regions to estimate 
the backgrounds both in section 111. A and lll.B, 
and take the difference as a source of systematic 
error. The result is 10.0% for the measurement of 
branching fraction and 4.2% for the upper hmit. 

4.8. Number of photons 

To estimate the systematic error from the re- 
quirement of two good photons, we compare the 
efficiency difference for this requirement between 
data and MC sample, and obtain 4.4%, which is 
taken as the systematic error from the two photon 
requirement. 

4.9. K* simulation 

The K* is simulated with a P-wave relativistic 

Breit-Wigner function BW — -, ^ifj J^", ^'i — 9, 

with the width r(s) = rpHa i+''V rP.i3 ^ 

is the interaction radius and the value (3.4±0.6± 
0.3){GeV/c)~^ measured by a K~n~^ scattering 
experiment [20] is used. Varying the value of r by 
Icr, the difference of the detection efficiencies for 
J/^P -> riK*°'K*°, J/t/j 7]Y{2U5) r]K*°'K*^ 
is taken as the systematic error from the uncer- 
tainty of the r value. 

The systematic errors from the different sources 
and the total systematic errors are shown in Table 
I. 

5. Summary 

With 58M BESII J/-0 events, the branching 
fraction of J/iJj riK*^K*^ is measured for the 
first time: 

Br( J/V) ^ TiK'^'K*") = (1.15 ± 0.13 ± 0.22) x 10"^. 

No obvious enhancement near 2.175 GeV/c^ 
in the invariant mass spectrum of K*^K*'^ is 
observed. The upper limit on Br{J/tp 
7yF(2175)) • Br(r(2175) ^ K*°'K*°) at the 90% 
C.L. is 2.52 X lO"**. Due to the limited statistics, 
we can not distinguish whether the y(2175) is a 
hybrid or quarkonium state. 
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